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Abstract. The singularities of the electromagnetic field are derived to include all 
the point-like multipoles representing an electric charge and current distribution. We 
show that for higher orders, it is more efficient to have fields represented in terms 
of symmetric and trace free moments. In the static case, the delta-singularities are 
expressed for arbitrary multipole orders, while in the dynamic case we restrict ourselves 
to the lower orders. The algorithm we give can be easily extended to the next orders. 



1. Introduction 



In the cases of electrostatic and magnetostatic fields of point-like dipoles, one has 
the well-known procedure of introducing Dirac (5-function terms for obtaining correct 
expressions of the electric and magnetic fields defined on the entire space. The 
corresponding field expressions take the following form [1]: 

1 , , 1 ( '^{u ■ p)v — p\ 1 



where ly = r/r, and 

-.^ / N 2/io ^, , «o f 3(1/ ■ m)iy — m\ 2iin ^, , 

B^{r) = J^m6{r) + g (^^ j = J^rn6{r) + (B)^.^, . (2) 

In these equations, by (. . .)r^o we understand an expression in which the derivatives 
are calculated supposing r 7^ 0, representing some well-known expressions of the fields. 
The expressions from equations ([1]) and ([2]) are introduced in Ref. [1] as conditions of 
compatibility with the average value of the electric or magnetic field over a spherical 
domain containing all the charges or currents inside. Another procedure for introducing 
equations ([I]) and (I2l) is based on an extension of the derivative didj/{l/r) to the entire 
space [2]: 

8<s4-(^^^)^^„-f%^W (3) 
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and this procedure will be applied in the present paper. A more pedagogical and suitable 
approach for understanding the origin of the difference between the electric and magnetic 
cases is done in Ref. [3]. Refs. [1] and [5] contain generalizations of the equations ([I]) 
and (|2]) to the dynamic case for oscillating electric and magnetic dipoles. 
Section |2] is dedicated to the reader less used with the tensorial formalism in handling 
multipolar expansions. For the informed reader, the section can be seen as an 
introduction to the notation and formulas used along the article. In Sections [3] and 
mthe results for the electric and magnetic fields are presented in the static case. The 
dynamic case is treated in Section [5l and the conclusions are outlined in Section |6l 
The formalism employed in this paper is a purely algebraic one. With a good 
understanding of the definitions and notation presented in Section [21 we think the 
reader will be able to verify every detail of the calculation and, possibly, to search and 
find some more adequate versions. 



2. Preliminaries 



2.1. Definitions, notation and formulas from the tensorial formalism 

Let us consider a n-th order Cartesian tensor denoted by T*^"-' and characterized by the 
components Tj^...j^, iq = 1,2,3. These components are, in fact, the components of a 
vector in the n tensorial product of the Euclidean space M^: 

^ = 'Til. ..in Gil ® ... ® G-i^ = Tjj...i,^ Gii...i„ ! 

where are the unit vectors of a Cartesian basis in M.^ and 6,^...,^ are the unit vectors 
of the tensorial product space. Some vectors from the tensorial product of spaces can 
be represented by the tensorial products of n vectors from the different space factors 
M.^. Particularly, in the case of identical factors, we employ the usual notation 

= ciii . . . cii^eii,,,i^ . 

The vector a can be the differential operator ^ = Cidi and, in this case, 

(V")=e.,...„9.,...9.„ . 

Because of this simplified notation, for avoiding confusions, the Laplace operator in the 

space M.^ is symbolized by A, and not by the frequent V^. 

We employ the following notation for the tensorial contractions: 

{^ii - i„_„iii---jm-^ii---im , n > m 
,n = m . 

A fully symmetric tensor S*^"\ which here is called simply "symmetric", has a projection 
on the subspace of symmetric and trace free (STF) tensors. Up to a numerical factor, 
this projection will be represented by the tensor 7"(S''"'*) = .S'-'"-'. For n = 2, for example, 
we can write 

S,, =5,,+<5.,A(S('') . (4) 
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The condition Sa = gives 

A = ^S,, . (5) 
In the case n = 3: 

Sijfc = <Sijk + 6{ijAk}{S^^^) , 

where by the symbol Ajj^ . j^j is denoted the sum over all distinct permutations of the 
indexes. The condition Suk = for k = 1,2,3 implies 

Though, maybe, only for a theoretical interest, let us consider the general case for 
the STF projection of the symmetric tensor S^"-* defined up to a numerical factor by 
the equation 

/\("-2) ^ A(S''"^) is a symmetric tensor and is defined by the condition that 7"^") 
is a trace-free tensor. For low values of n, the ones of really practical interest, the 
components Aj^ . j^ ^ can be calculated directly from the equation system representing 
the vanishing relations of all the partial traces of the tensor 7"'^"^ However, we mention 
here a general formula known from literature [6l [7] which, with the notation from the 
present paper, is written as 

n-rcWii _'^' (-ir(2n-l-2m)!! ^ 

m=0 ^ 

The symbol [a] represents the integer part of a and denotes the components 

of the (n — 2m)-th order tensor obtained from S*^""-* by contracting m pairs of symbols i. 
This equation is known as the detracer theorem [7]. As a consequence of this theorem, 
the components of the tensor A^'"~^^ are written as 

A rc(n)i_'"v^" (-ir[2n-l-2(m+l)]!! 

'Hi...J„_2[-3 J 2^ (m+ l)(2n - IV^ "{*1«2 ■ ■ ■ "*2m-li2m^.t2m+l...in-2}- 

■m=0 ^ 

These formulas are useful for defining and processing the multipole expansions of the 
electrodynamic potentials and fields. 

2.2. Multipole expansion of the electromagnetic field in Cartesian coordinates 

The multipole expansions of the potentials in the Lorenz gauge are written as 

<|.(r,^) = -^V^^V"||^-^ (6) 



and 



4:71 ^ n\ 

n>l 



VxlV'-'II^^Vv-'ll?^ 



(7) 
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The dot symbolizes the time derivative and r = t — r/c is the retarded time with respect 
to the origin O of the Cartesian axes in the point corresponding to the vector r . The 
origin O is a point from the support of the electric charge and current distribution. The 
tensors P(")(t) and M(")(t) are the electric and magnetic moments of the electric charges 
and currents distributions p{r,t) and J{r,t), the support of these distributions being 
included in the domain V: 



p(")(t) = f d3xr"p(r,t), 
Jv 



and 



= I r" X J{r,t) . (8) 



n + 1 

In the last equation, a tensorial contraction via the Levi-Civita pseudo-tensor Sijk is 
employed: 

which in the particular case of T*^") = becomes 

X a = . . . b,„_, {b X a) e,,...i„ . 

The expansions ([6]) and ([7]) are running in the exterior of the minimal radius sphere 
including the support of p and J. 

From equations ([6]) and ([7j) one obtains the following expansions of the fields E{r, t) 
and B{r, t): 

^/ X ^ dA(r,t) 1 V- (-1)"-^ P^^Hr) 
E(r,t)= -V$r,t = 1 > ^ — — 



/iO (-1) 



47r ^ n! 

n>l 



(it 47r£:o 

n>0 

Vx|V»-'II^Vv»-'ll^!^ 



(9) 



and 

S(r,t) = V X A(r,t) 
/io ^ (-1)"-^ 



47r n! 

n>l 



<10) 



For the magnetic moments ([8]), it is also possible to introduce STF moments 
_^(n) _ t'j^mC")) this time, there are two steps required in order to complete 
the objective. The tensor M*^*^) is symmetric only in the first n — 1 indexes and satisfies 
the property 

Mii...j„_2 gg = . 

In the first step we must obtain the symmetric projection (up to a numerical factor) 

M of the tensor M We begin with the first simple example corresponding to n = 2. 
Let us write the identity 

1 1 -^1 
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^ (2) 

where M is the symmetric part of M*-^^ and 



N,(M(2)) = e,,uM,k = 11 d=^a; [r X (r X J)] 

<J Jv 



In this case (n = 2), Al'^-* =M and, consequently, corresponds to the STF projection. 
Therefore, 

M., =M, + ^£^ifcNfc(M(2)) . (11) 
For n > 3, we can generahze this resuh writing the identity: 

n-1 



n 

A=l 



where by nI^^"* . ) we understand the component without the index ix and the tensor 
|\j(n-i) ^ m(mW) is given by 



(M^ ^) — Si^_^pq Mi^...i„-2Pq ~ I 

It is a tensor of the same type as M'^"~^\ i.e. symmetric in the first n — 2 indices and 
with n — \ vanishing traces {^i^...i„_3qq = 0). Therefore, the STF moment Al*-"^ is given 
by the components 

M^,,„i^ =M,,„,i„ -<5|,,,2^3...i„}(M ) . 

As it will be seen in the dynamic case, we have to express even the symmetric projection 
of the tensor N(M^"'^). For this, it is useful to introduce the operator M defining the 
correspondence 

Repeating this operation, we obtain: 



AA2/c^M(")) = ^il^ I d 

n + l Jv 



3^ ^2k ^n-2fc ^ J 



Jspk+ln^{n)^ ^ ( ^) ^ /" r^'^ x (r x J) . 

^.c?. Delta- function identities and multipole singularities 

Usually, the multipole expansions are considered, term by term, as functions of r defined 
on excepting a singular point which is chosen as the origin O of the Cartesian 
axes. Actually, these multipole terms are mathematical distributions (or generalized 
functions) considered firstly as regular ones having as support the entire space except 
the origin point O. The observable quantities are expressed as weighted averages on 
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spatial regions or as surface integrals of functions of field variables. No problems appear 
when in these regions r 7^ 0, but when we have to calculate for example the interaction 
of a distribution (p, J) with an external field {E^ B), 

W,nt = J d'x (p$- J- A) 

and the system associated with the external field {E, B) is represented by a point- 
like multipole system placed in O, the singularities of the potentials or fields having as 
support this point become unavoidable. As done in [2], these singularities are determined 
starting from some ^-function identities associated with the extension of multiple spatial 
derivatives of the functions of the type 1/r to the entire space. Having in mind the 
dynamic case, too, we generalize such identities to the derivatives of the function fir) jr. 
Representing the corresponding derivatives as functions of r and t for r 7^ 0, we can 
write their expressions in the form 

1=0 

where C\^'^]^ are fully symmetric in the indexes The general form of the 

coefficients C is a simple consequence of the symmetry properties and of the derivative 
rules: 

[fl 

'^tX = Yl ^^"''^ "^{^1*2 • • • Si,,^,i2k ^2k+l ■■■Vn} ■ (14) 
fc=0 

For lower n the coefficients C and K can be easily calculated by successive derivative 
operations. As one can see in the following, for obtaining general results for the 5-type 
singularities of the multipole terms, it will be necessary to know only the coefficients 

= (-l)"(2r2- 1)!! . (15) 

y(nO 



In Appendix A we enumerate the coefficients for the first five values of n. 

Considering the distribution representing the derivative ( fT3l) extended to the entire 
space, we understand by (9j^ ■ ■ - dt^ if {'^) / ''^)) (o) the singular part having as support the 
point O. In the following, we are interested in writing explicitly this contribution. 
A correct procedure for extracting this singular part is that employed in [2] defining 

fir) 




(Txdi^ . . . di. 



r 



(16) 



In this equation, 0(r) is supposed an element of the domain of distributions i.e. an 
infinitely differentiable function. As it can be seen from equation f lT^ , similar restrictive 
properties are considered for /(r). The domain is the spherical region delimited by 
the spherical surface with radius e. A basic hypothesis is the existence of the limits 
considered here, but the general demonstration of this property is a purely mathematical 



On the 6 -singularities of the electromagnetic field 7 

problem and bypasses the objectives of the present paper. Therefore, the result of the 
limit considered is independent of the domains defined only by the condition — )■ O. 
Finally, for the calculation process in the present paper, the essential property is the 
possibility to represent the functions / and in the domain by their Taylor series 
upon the origin O: 



A=0 ■ ^ / r=0 

Equation (fT6l) can be written as 

( (Ai..i„(/))(o) , 4>) = Imi i dSi^i.di, . --di,^^ - ( (A2..i„(/))(o) . ^n^) 
generating a recursive calculation for a given n. 

Let us express the (5-singularities for lower order derivatives. Even in these simple cases, 
it becomes obvious the importance of the angular average of a function g{u) defined as 

For this average, we have the well known formula [6]: 

z/, ) = { 1 °' n = 2k + l, ^^g^ 

[ (n+l)!! ^{hh ■ ■ ■ ^in~iin}^ U = 2k, /c = 0, 1, . . . 

The first 5-singularity is obtaining for n = 2 calculating the limit 



lim 



The first surface integral limit can be written as 

L,, = hm / dn{iy) V, [-^/(r,) - ^/(r,)! [0(0) + ev^ {dk4>\ + ...], 

where = t — e / c. Since /(rg) = /(t) + 0{e), and all the terms containing as factors 
positive powers of e vanish at the limit e — )■ 0, we can write 

Xi. = -47r {u,u,) f{t)m = 5,, /(t)0(O) . 

The second surface integral limit 

= lim ^ dn{iy) evj f{re) [(<9i0)o + euk {didkct))^ + ...]= , 

since all the terms contain as factors positive powers of e. As the volume integral limit 
also vanishes, we can finally write: 

=-^fm,5{r). (19) 

^ / (0) 
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Including the expression of the derivative for r 7^ 0, the above result written for f{t) = l 
verifies equation (j3]). 

Let us consider the next order derivative and calculate 



(A,.(/))(o) , 0) 

lim |y"dfi(z/) 



= lim 

4^ ' .3. 



X (0(0) + eui (9,0)o + ...). + —/ d^x f{t)6,kS{r) d,<P{r] 

where equation ( |T9l) is employed in the volume integral. In the surface integral limit, as 
e — )■ 0, the terms containing as factor the second order time derivative / vanish because 
of the positive powers of e. The terms containing the first order time derivative / are 
proportional to positive powers of e for the term corresponding to 0(0), but in this 
case the corresponding limit cancels since it contains also as factors the null averages 
{yiVjVk) and {vi). The only terms giving a non vanishing limit are provided by the 
product /(t^) {dicf))^ such that we can write 

(Ajfc(/))(o) ,0/ = 47r {SuiUjUki^i - 5jkViVi) f{t) {dicj))^ + —f{t)8jk {d,(f))^ 

(\ \ \ 4:71 

= 4n I -S{ijSki} - -SjkSiij fit) {di4>)^ + Y/(t)5jfc (9,0)o . 



Finally: 



d,d,dJ-^\ = -^f{t)6{,,d,}6{r) . (20) 

/ (0) ^ 

For f(t) = 1, equation ( l20l) becomes the 5-singularity corresponding to equation (13) 
from Ref. [2]. Equations ( IT9|) and ( 120|) will be applied in the static case for f(t) = 1. 

3. Singularities of the electrostatic field 

Let us consider the multipole expansion of the electrostatic field derived from the 
potential expansion ([6]) in the static case: 

E(r) = tll!l^V"+H|^ = izll!lipW||V"+ii 

47r£o^ n! r 47r£o ^ n\ r 

In the dipolar case {n = 1), the potential has no 5- singularities. For the electric field, 
we apply equation f lT9|) obtaining the known result ([1]). 

A first interesting feature of the (5-singularities problem appears beginning with the 
electric quadrupole potential and field. The singularity of the potential is obtained 
employing equation ( !T9|) : 

/I P(2) \ 1 

VSvreo r J Geq 
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The first term from tlie riglit-liand side of this equation is invariant to the substitution 
p{2) _!. -pCs) ^ 'y-(-p(2)^ ^^^^ employing equations @ and ([5]) with S = P^^\ we can 
write 



Usually, since in the exterior of the charge distribution the potential $(r) is solution of 
the Laplace equation, it is represented by a series of spherical functions: 

oo ^ I 

^(^) = E ^ E ^) • (21) 



«=0 ^ m=~l 



The 2/ + 1 spherical moments corresponding to a given / are linear combinations of 
the 2/ + 1 components of the l-th order STF tensor 7"-'^ For the above example of 
the 4-polar potential, it is obvious that using the expansion representing the particular 
case of the equation f l2T]) . the cie/ta- singularity of the potential is lost. The conclusion 
can be extended to the higher-order multipolar potentials. For the electrostatic field 
E{r), the employment of the corresponding expansion in spherical functions leads to an 
incomplete description of the (5-form singularities. 

Based on this observation, in our opinion, some care is necessary when passing to the 
representations by spherical functions since it is equivalent to working only in a subspace 
of the mathematical objects used correctly in the theory. The direct substitution 
P — or passing to the spherical tensorial representations is not always justified. 
This observation is pointed out in another context in [12], too. 

Searching the ^-singularities of E^'^\r), instead of a straightforward processing of the 
corresponding expression in terms of the primitive moment P*^^-*, we can introduce from 
the beginning the STF moment 7-*^^^ writing 

1 / P(2)\ 1 / 7?(2) A 



(°) Sirso V ^ y (0) ^^^0 V ^ / (0) 

1 / 7^(2) \ 1 



Sireo V ^ / (0) 2£o 
since A(l/r) = —A7r6{r). Equation (!20ll finally gives 

E^'\r) = {E^'\r)),^ + J-p(2)||V5(r) + -LAV5(r), A = A(p(2)) , (22) 

where 

1 /_.,,P(')A 1 



Even from this simple example it becomes obvious that this calculation version is 
convenient for higher n since some contractions of T'^") with the Kronecker symbols 
vanish. 

In case n = 3, the advantage of inserting from the very beginning the STF moments 
is more evident. Instead of searching the singularities of the fourth-order derivative 
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didjdkdi (1/t), we search the singularities of the contraction of this derivative with the 
STF moment 7"^'^^ writing 

77(3)- 



Cj hm 

(0) / 



For the surface integral limit, we apply equation f lT3|) particularized to the static case. 
Denoting by L^r this limit, 

= e, lim £ ^^^^z/, CffJi^^Vjki 0(r) 
and, inserting equation (]A.4p . one obtains 

= Bi hm ® — [-IbuiUjUki^i + Sz^idijTct'/}) Vjm 4>{r) . 

Since S{jki^i}Vjki = 0, we get 
La= - ISci lim / UiyjVk^iVjki (0(0) + ez/q (9g0)g + ^eV^z^, (5g5.0)o 



- 47r X 15 lini ^ ^UiUjUkUi Vjki (^0(0) + ez^g (<9q0)o + ^ 



Denoting by a the power of e in the Taylor series of 0(r), we easily see that for a > 3 
all the terms from the last equation vanish for e — )■ 0. For a = 0, the corresponding 
term vanishes since {i'ii'ji^k^i)Vjki = 0, as it is seen by inserting equation ( IT8|) . For 
a = 1, the related term contains the average of an odd number of factors v yielding a 
null result, too. The only non vanishing result corresponds to a = 2 such that 

15 

= -An X Y Si (uiUji^ki^iUqiys) Vjki ■ 

In the last equation we have the contraction {i'i^i'i2i'i.^i'i^Vi^Vi(^)Vi^i^_^i^. It is not necessary 
to write explicitly the angular average of the six factors v since, considering the general 
form of this expression given in equation (fT8|) . we see that only the 3! = 6 terms of the 
form 5i^i.5i^i^,5i^i^ with the indexes j, A;, / covering all the permutations of (4,5,6) give 
non vanishing limits. Consequently, 

= (a,a,0)„ = (VV)o ■ (23) 

The volume integral limit is zero since 

h— =0, 
^ Ao) 

as it can be easily verified. Equation ( 123|) implies 

24 
14 



7^(2) \ 24 

V^ll^ =_flp(3)||v25(r.) . (24) 



(0) 



In [Appendix B| this equality is proven for the general case. Let us write 

(E^'\r)),^, = ^— ( V^W^] = Lp(3)||v2<5(r.) - J-A\\v'S{r) 

^ ^ ^no) 247r5o V r ) lAe^ " ^ ^ 2£o " ^ ' 
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or 



where 



('■)), 



1 .v^iifi!! 



1 -v^ll^ 



Even if only for the mathematical interest, we mention the possibility of expressing the 



multipole electrostatic field for arbitrary n. In [Appendix B| it is introduced the general 
relation corresponding to equation f l24|) : 



Ann 



(0) 



2n 



pW||V"-i5(r) 



such that 

(^;(")(r)) 



\n-l 



(0) 



A-keq n\ 
(-l)"n 



7n+l 



.-pW _|_ 



n 



rn — 1 



(0) 

5(r) 



(25) 



£o V^'(2^ + l) 2ra! 
We point out that for n > 4, for obtaining the complete singular structure of the 
field, we have to introduce the irreducible tensors in the expression from equation 
corresponding to the (5-singularities . This implies a recursive calculation introducing 
in equation ( !25|) the STF projections of l\f^~'^\ JK^^~^\ .... For example, in the case 
n = 4, we have to reduce the tensor A*^^) = A(P'^^)) writing: 

A(2)||v3(5(r) = r(A(2))||v3(5(r) + ^{^^^^) VA6{r) , 

where 

1 



{2)^ 



A(A 



7-,(A(2): 



-P 

20 ' ggss ; 



^ P 



-P • 

g' qqss'^i] 



4. Singularities of the magnetostatic field 

Let us write the multipole expansions of the vector potential A{r) and of the magnetic 
field Bir) in the static case: 



and, respectively, 

B{r) = ^y 



n>l 



n\ 



-V X (V""^||M(")) 



/io (-1 



n>l 



V V 



A V 



(26) 



For r 7^ 0, the multipole magnetic field B{r) is expressed, as in the electric case, 
only by a scalar potential since the second term from the right-hand side of equation 
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vanishes. However, the 5- singularity of the magnetic field is not similar to the 
corresponding singularity of the electric field, this second term giving the difference. 
The difference becomes easily obvious if one considers the fictitious magnetic shells (or 
sheets) employed in the Ampere formalism. It suffices to consider the case of the point- 
like magnetic dipole which can be taken as the limit of a current loop of infinitesimal 
size. For finite dimensions, the field of the loop is derived from a scalar potential $m 
which is defined by an integral on the corresponding sheet and having a "jump" in all 
their points. Just this jump generates the (5-singularity corresponding to the second term 
from equation ( 126|) . An explicit calculation of this limit when the loop concentrates in 
a point is given in Ref. [T3] . 

In the dipolar case (n = 1), retaining the singular term, we have 



(0) 

Processing the first term as in the electrostatic case, we obtain the well-known expression 
(I2]). No (5-singularity is present in the dipolar vector potential. 

Let us consider the quadrupole term from equation (l26l) and the corresponding 
singularities: 



(0) 



« (2) 

By introducing the STF moment M.^"^' =M defined by equation ( fTTI) . 

= _ ^A4(2)||v5(r) -^Nx V5{r) , 
10 II ^ ^ 4 ^ ' ' 



(0) 



where N = e,;N,;. It follows that 

^'r^o 10 II ' ^ 4 



5(2) (r) = (B(2)(r))^^^ - ^ M(2)||V5(r) - ^ iV x V5(r) (27) 



where 



Let us consider the quadrupolar vector potential A^'^\r) with equation ( 1TT|1 already 
inserted: 

A<=)(r) = - ^e,e,,a,a^ = V x f V||^ 

8tt r Stt \ r 

lOTT \ r / lOTT r 

The term proportional to A(l/r) vanishes for r ^ 0. Whereas B^'^^ is invariant to the 
substitution P^^^ — )■ this property is not verified by the vector potential A^'^\ 

This substitution has as result an additional gradient term, i.e. a gauge transformation 
of the potential. Though, by the extension to the entire space, the gradient term 
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from expression generates a 5- singularity, this term gives no contribution to the 
singularities of the magnetic field The singularities of B are, by definition, only the 
singularities corresponding to the curl of the vector potential. It remains still the 
problem of independence of the physical results on the A (5-singularities corresponding 
to different gauges. Maybe, the presence of such singularities in the expression of an 
interaction Hamiltonian corresponding to the density term J ■ A can generate the 
problem of proving the gauge independence of physical results as in the case of the 
Aharonov-Bohm effect. 

We can give, as in the electrostatic case, the general formula for the singularities of the 
magnetic field writing 

(flC,;,) , fv....||^ _ V-||a!^) . (29) 

^ ' ^(0) 47rn! V ^ ^ J (q) 

The introduction of the symmetric moment M , equation ( |T2l) . gives 



M(") M 1 K ^ M 



, , t\tnq " I I 

n — ' r r 

A=l 



A=l 

since in all terms of the sum there are the contractions s^i^qd^di^ for / = 1, ... — 1 
which vanish. From equation ( l25l) we can write a corresponding result for the magnetic 

field with the substitutions P(") , P^") A<("), Eq l//io: 




_A4W||V"^i(5(r-) 



where A = A(M ). The second term from equation fl2^ is processed as follows: 

'VI „ 1,, . IVI 1 „ „ ■^-^ . (li...ln^iq) 



n ^ — ' r 

A=l 



O (n) 



= V '\\A + eieiq,^_, N,,...i_2A • ■■dir.-A- ■ 

r n r 

■H- 

The term containing M can be further transformed to: 

M^"^ ||V"-i5(r) = M(")||V""i5(r) + e, 9,, . . . 5(r)5|,,,,A,3„,„_,,| 

= A4(")||V"-i5(r) + ^^^^l|^^e,9, . . . ^,^_.AKr)K..^.-^^ 
+ {n- l)eidi di,... di,^_^5{r)\,„i^_^ 
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Finally, 

^ ^ ^^(0) n! [V 2n + 1 2 y " ^ ^ 

where the notation 

(a*- I X IB*- ^) j^^^ = £i„qs^ii...i„-.i q^ii...in-i s 

is introduced. The general expression is written as 

SW(r) = (BW(r))^,^„+(sW(r))(„) (r) , 

where 

(B(«)(r)) = (-l)"-Vo vn+i||-^^"^ = (-irVo ^„+,,,M(") ^ 
'^^■t^o 47r?T,! r 47m! r ' 

with the derivatives calculated for r 7^ 0. 

Here is the place for the same observation as at the end of Section ([3]), but in the 
singular term containing N*^""^-* the reduction to the STF tensor T*^"~^) begins from 
n = 3. Particularly, 

<->(2) <_> 1 



and 



r(N(2)) = N , N - (N,, + N,,) , 

<->(2) 1 1 

Obviously, the analysis of the higher order singular terms is more complicated for the 
magnetic field than for the electric one, but it is purely a technical problem. 

5. The dynamic case 

Considering the n-th order multipole electric and magnetic fields, writing separately the 
contributions of the electric and magnetic multipoles, 

' 4:1160 n\ \ r r 

E(-'-^)(r,t) = , V X ( v"-!!! *^ ^^A (31) 



and 



BK^)(,.,) ^ (-1)"-;^. fv....||M^ , V'-||AM!!!(li ) , (32) 



47rn! 
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we search the corresponding 5- singularities. Processing the different terms from equation 
f l3T]) and f p2|) for reducing the moment tensors to the STF ones, we deal with the Laplace 
operator A applied to functions of the type /(r)/r. Since the equation verified by these 
functions is 

/(r) 1 f{r) 



A- 



-47rf{t)6{r) 



(33) 



r dt"^ r 

the processing is more complicated in the dynamic case compared with the static one 
For the electric field of the electric dipole we write the singular part as 



(^;(i'P)(r,t)) 



(0) AttEc 



V V ■ 



(? r 



(0) 



Atteq 



did. 



(0) 



Note that the second term proportional to p has no (5-singularity. Equation f[T9|l gives 
a result similar to that from the static case: 

(i;(i'^)(r,t))(„) = -^p(t)5(r). 



3£:n 



The magnetic field B^^'P\r,t) of the electric dipole has no 5-singularity. 

Obviously, the electric field E^^'^\r^ t) has no 5-singularity either. Writing the singular 

part of S(i'™)(r,t), 

1 rhir) 



(B"'""('-.')),o, = £ 



V V- 



^2 y. 



+ A'Km{t)6{r) 



(0) 



one obtains the result similar to that from the static case: 



(B(i'-)(r,t)), 



2/io 



m{t) 5{r) 



'(0) 3 

Let us consider the electric field E^'^'P\r,t) and search the corresponding (5-singularities 
by introducing firstly the STF moment T"^^^: 

P(^)(r) 

"^2 



(i;(2'P)(r,t)) 



(0) SttEo 
1 



(0) 



V A 



QTlEo 

The insertion of equation 



p(2,p) 



A(r) 



c r 



(0) 



gives 



Stieq 



-47rA(t)V5(r) 



(0) 



Expressing the singularity of the first term from the parenthesis by the insertion of 
equation (l20l) . we can finally write 



^;(2'^)(r,t)= (£;(2.^)(r,t))^^^ + ^p(2)(t)||V5(r) + ^A(t)5(r) , (34) 



where 



(£;(2'^)(r,t)) 



Stt^o 



ivii^ 
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The result fl34l) is similar to the result f l22p from the static case. For the electric 4-polar 
term, the electric field expression for r 7^ is invariant to the substitution P^^^ — 7^^^^ 
as in the static case but, as we will see in the following, such property is not yet verified 
for higher order multipoles in the dynamical case. 

Searching the singularities corresponding to the electric field EJ^^'™) of the magnetic 
quadrupole, the difference from the static case is obvious. Let us express this field with 

the help of the STF magnetic moment Al^^^ 



«(2) 

= M : 



i;(2'™)(r,t) 



+ -e-iSijkSikqdjdi 



N,(r) 



V X VI 



M 



(2)- 



1 



V X 




We employed the same notation as in equation (p7|) . The singularities added by the 
extension of this expression to the entire space are the following: 



V X VI 



(2)- 



.Nir) 



(0) 



(0) 



djdi- 



Mi, 



J (0) 



An 

■ — eieijkMjk{t)Kf) = 



-^N{t)S{r). 



(35) 



For the last term from the expression of E^'^'^\ we have to consider the equation (133|) . 
i.e. 

^Mt) 1 N ir) 



r c 
such that we can write 



i;(2'-)(r,t) 




1 N (r) 

4^4 ^ 



+ 



;N{t)6{r) 



1 



;N{t)5{r) 



12£oc2^''^'"''' Ae^c^^'^"'"'-' ' ' ^^^^ 
In the above result one can notice that for r 7^ 0, the expression of the electric field is 
not invariant to the substitution P*^^^ — 'P^'^\ The additional term which appears 
as a consequence of this substitution is exactly the expression of the electric field 
corresponding to an electric dipole with the moment 



(3T) 



In equation (1361) . we have written separately the term which represents the 5-singularity 
corresponding to this dipole. 

Let us express the field B^'^\r,t). For the part B^'^'P\ we insert the STF moment 7-*^^^ 
which gives 




.(2), 
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(3^ 

since SijkdjdiSik = Sijkdjdk = 0. It follows that S^^'P) is invariant to the substitution 
p{2) _^ j>i2) g^j^^ j^g^g ^-singularities due to equation fl35l) written for Al^^^ — T'^^^ 

<^(2) 

The introduction of the STF magnetic moment Al*-^-* =M in the expression of B^'^'^"'^ 
results in 

B(.»..(,,i)=_^^lfv^.||M!!M_v||A^"'M 



Stt \ r r 

— ^ - V A — + -V X A- ^ ' 



Itx \ r r 2 r 

Separating the ^-singularities with the help of equations ( 120|) and (133|) . we obtain 

,(2), 



Stt \ r c"' r 2c^ r 



- ^M(^)(t)||V5(r) - ^iV(t) X V<5(r) . (39) 
Since for r 7^ the magnetic field of the magnetic quadrupole can be written as 

we can conclude that for r 7^ this field is not invariant to the substitution M^^^ — )■ 
Al*-^^ and the additional term introduced by this substitution is equivalent to that 
of an electric dipole having the moment 5p' given by equation ( 137|) . Obviously, the 
corresponding term has no (5- singularity such that it is not represented by a singular 
term in equation (l39ll . 

Let us go further to the field of the third order electric multipole. Beginning with the 
electric field, we have to extend to the entire space the expression 

B<-.(,,() . J- f V'll?^ - iv^ll^) . (40) 

^ ' 247r5o y r " r ^ ^ ' 

We consider separately each term from the last equation introducing the STF moment 
p(3). 



ty ' ' 

Inserting equation ( 1B.8I) for the (5-singularity of the first term from the right-hand side 
of the last equation and employing equation (l33l) for the second one, we can write 



2 



4-7r 

-V^'-\t)\\VH{r) - —k{t)5{r) - l2nK{t)\\VH{r) , (41) 



C2 
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where A = AjCj. Regarding the second term from equation (HU]) . we similarly obtain 



V^||P!!M_(v.||5!V)^2V^||A(l)+iMl) _?2£Amr) (42) 
rl r r c r f 3 

since the first term containing V^||('P^^^(r)/r) has no (5-singularity. The insertion of 
equations (HTl) and ( H2ll in equation ( HOl) gives 




14eo 2^0 9eo 

From the above equation we can conclude that E^^'P\r,t) is not invariant to the 
substitution 7-*^^^ — )■ 7"-^^ This substitution introduces an additional term which has 
precisely the expression of the electric field of an electric dipole of moment 

Sp" = 7^A(t) . (43) 
Let us express the magnetic field B^^'^^: 

S(3'^)(r,t) = X ( V'W^^] = -^e^e.^.dAA.^^'^ -(44) 

V ' ^ 247r \^ " r J 247r ' ' r ^ ^ 

The introduction of the STF moment 7"-^-* has as result 

^ ' ^ 247r y " r ^ 247r ' ' r 

and since V x (V^IK'P {t)/^)) has no ^-singularities, we can write after simple 



algebraic calculation and after using equation ( 1331) that 

,(3) 



^ ' ^ ' 247r ^ " r 47r Gc^ r 



0A(t) X V5(r) . 



For r 7^ 0, the additional term introduced by the substitution P*^^^ — )■ T-*^'^-' in equation 
(jSj) represents the magnetic field corresponding to the electric dipole with the moment 
dfl defined by equation P3l) . 

Let us consider the electric moment up to n = 3, and the magnetic one up to n = 2 
and the corresponding sums of fields representing one of the first approximations of 
a complex system assimilated with a point-like multipolar system. From the above 
results, as it was expected, we see that the substitutions of the primitive moments by 
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the corresponding STF ones introduce some additional terms. The sum of these terms 
is equivalent to the substitution 

P -> p = p + Sp , 

where 

6p = 6p' + Sp" = -\i, 
and the vector t is defined as 

t = In - -A = — [ d^x iir.J) r - j) . 

This vector is the so-called electric dipolar toroidal moment and represents a first term 
from a series of toroidal moments introduced by Dubovik et al [13]-[T6] by generalizing 
Zeldovich's original idea that a closed toroidal current represents a certain new kind of 
dipole [T7] . The corresponding (5-singularities can be easily identified in the expressions 
established above. 

The procedure can be continued for the next orders, the technique being the same as in 
the case of the first orders treated above. Unlike in the static case, in the dynamic one 
it is more complicated to establish expressions of the singularities for an arbitrary n. 
The recursive character involving a general number of n steps makes such relations hard 
to derive and apply [11]. Furthermore, at the current level of applications, we believe 
that it is mandatory to have expressions up to n = 3. 

6. Conclusion 

In this paper we have expressed the (5-singularities of the electric and magnetic 
multipoles. The results were given for arbitrary multipole orders n in the static case, 
while in the dynamic one, the outcomes correspond to the lower orders. However, the 
calculation algorithm is presented in a manner that facilitates the processing of the 
results step by step to the next orders of the dynamic expressions. The central idea of 
the article is that, instead of employing the ^-singularities of the functions /(r)/r, it is 
more efficient for the higher orders to search directly such singularities for the multipole 
fields represented in terms of the STF moments. Working in Cartesian coordinates and 
employing a particular system of parameters and notations adapted to the technique 
of reducing the primitive Cartesian moments to the corresponding STF ones, we have 
stressed the significance of the first type of moments in the process of searching the field 
singularities. 

In our opinion, a first field of application of these results is to extend the classical 
argumentation of the generalization of the Fermi contact term to arbitrary multipole- 
multipole interactions, useful in the studies of the hyperfine interaction, [HI [T9] . 
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B B /(^)-V ^ ^(n.0 d"-7(r) r 

C(°'°) = 1 . (A.l) 
Cf^'°^ = -.., Cp^ = -^. (A.2) 
= J^iJ^j , C'j^/ = Sz/jZ/j - Sij , C^^'^^ = SuiUj - Sij . (A.3) 

C(3'2) = _ I5z/,z/jz/fc + 3 6{ijUk} , C^^'^) = - IhuiUjUk + 3 . (A.4) 

^ijid = 45 ViViVk^i - QS{ijJ^kJ^i} + %4/} , 
Cf/fc? = UiVjUk^i - 15 + 3 , 

Appendix B. Some general formulas for STF tensors singularities 

When searching the ^-singularities of the electromagnetic fields, we have to calculate 
singularities of some derivatives of STF electric and magnetic moments as 

^n+iu}__Xn^ Vx ( V"-i||-!-^), (B.l) 
r \ r J 

where T'-"'-* is one of the STF electric or magnetic moments. 

Let us search the ^-singularities of the first expression for n > 2 starting from the 
definition 



V"+i||-l ^ 1 ,(/))= lime,- 



(0) ' 



d\9,,...9,„I^^^0(r)9,0(r) 



(B.2) 



Denoting L^j the limit corresponding to the surface integral from the last equation, the 
insertion of equations f|T3|) . f|T7|) gives: 

1,.= Ihne, 



/=0 a>0 A>0 



C jn-Z+A 

i dfi(^.) z/, Ct::1j^^n ■ ■ ■ ^j.T^^jH...Ut) {d,, . . . B,^ 0)o . (B.3) 



X 
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Considering the expressions dUD of the coefficients C^"''), in equation f lRSjl we can 
perform the substitutions 

ll...ln U ''1 ''Ti 

without changing the result since the contractions of Tj^ with all the terms from 
the coefficients C\^'^}^ containing at least a symbol Kronecker vanish. Expressing the 
integral in this equation as angular average, we can write 

(_l)A^a-/+l+A^(" 



4:71 e,- lim 



■ "A"o ^ ^ ^ c"~'+^a;' A' 

/=0 Q>0 A>0 



n-l+X 



T,„„,„(t)(9,,...9,„0)o . (B.4) 

As we can see from equation (1181) . the contractions 

are different from zero only if a + 1 + n > 2n, i.e. 

a > n-1 . (B.5) 

Let e = a — / + 1 + Abe the power of e in equation (lB.4p . Since / < n and A > 0, the 
inequality fIB.SP implies e > 0, the equality holding only for a = n — 1, I = n, X = 0. 
-Lo- 7^ only in this case since for e > the limit vanishes because of the positive values 
of the powers of e. Therefore, equation f lB.4p becomes 



4n e,- 



^-^-^ {ui^... Vi Uj^ . . . Ti,...i„(t) [dj, . . . dj„_, 0)g . (B.6) 

As one can easily see from equation f[T5]l . the remaining contractions are different from 
zero only for the n\ terms of the form . . . (5i„_ij„_i5i„j, and, inserting the value (fTSjl 
of Kq"''"'-' we can finally write 

The volume integral limit from equation ( IB. 21) can be processed applying repeatedly the 
Gauss theorem. The new surface integral limit can be treated as in the previous case 
since it is represented by an equation similar to equation ( 1B.4P with n — )■ n — 1. The 
contractions to be considered this time are 

They are different from zero only if a > n and, since / < n — 1, the exponents of e are 
e = a — l + l + a > 2 such that the limit of this integral for e — )■ cancels. Therefore, 
equation (IB. 71) inserted in equation (IB. 21) has as final result 

/^„+1||T5^\ =4^ (~^)"^ T(")(t)||V"-^^(r) . (B.8) 

V ^ y (0) 2n + l 

The curl expression from equation (IB. 10 has no ^-singularities and this result is obtained 
applying the same procedure as in the case of the ffist expression from this equation. 
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Indeed, processing the corresponding surface integral limit, we can reduce it to terms 
having as factors the following contractions 



But 



therefore, 



^ijk {^ii ■ ■ ■ ^in-i ^ji ■ ■ ■ ^jn-'l) Tii...j„„i fc(^) ~ ^ijkJ i\...%n-ijk — 



V X V 



. 



(B.9) 



(0) 
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